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Section  1 


INTRODUCTION  AND  SUMMARY 


This  report  deals  with  SAI's  continuing  program 

to  rigorously  characterize  photoemission  for  soft  X-ray 

sources  incident  on  various  materials  found  on  space 

systems.  Over  the  past  three  years,  DNA  has  jointly 

sponsored  the  given  theoretical  effort1-10  and  an  experl- 

11-13 

mental  one  by  Bernstein  of  ARACOR  to  define  photo¬ 

emission  properties  for  plasma  radiator  sources  (PRS). 
Such  information  allows  SGEMP  analysts  to  model  the  elec¬ 
tromagnetic  environment  of  systems  exposed  to  PRS  radia¬ 
tion  in  simulation  tests.  Overall,  the  joint  program  has 
been  a  successful  one.  Theory  has  provided  a  detailed 
description  of  the  features  found  in  photoemission  as 
well  as  a  detailed  description  of  the  transport  process 
itself.  Generally  good  agreement  has  been  achieved  be¬ 
tween  theory  and  experiment  in  the  primary  photoemission 
yield  for  most  materials  considered.  Such  agreement 
provides  confidence  that  photoemission  is  being  well 
characterized.  Where  agreement  has  not  been  so  satis¬ 
factory,  further  investigation  has  been  carried  out 
leading  to  better  understanding  of  the  given  material. 

PRS ' s  do  not  provide  as  much  control  as  one 
would  like  to  have  for  comparing  theory  with  experiment. 
Not  only  are  there  variations  in  the  X-ray  output  from 
shot-to-shot  but  there  are  also  the  rapid  time  varia¬ 
tions  within  a  given  shot  on  the  scale  of  tens  of 
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nanoseconds.  Continuous  low  fluence  X-ray  line  sources, 
on  the  other  hand,  in  principle,  provide  for  complete  con¬ 
trol  of  experimental  conditions.  Based  on  these  consider¬ 
ations,  with  the  desire  to  further  test  photoemission 
characterization  from  the  joint  program,  DNA  has  sponsored 
photoemission  measurements  by  Bernstein  for  X-ray  line 
sources  and  predictions  by  us  using  code  SXRP  for  the 
same  sources.  Many  of  these  theoretical  results  will  be 
presented  in  this  report. 

DNA  has  been  the  primary  sponsor  for  development 
of  code  SXRP.  This  is  a  large  code  which  solves  the 
Boltzmann  transport  equation  in  matrix  form  derived 
from  the  integro-diff erential  form  of  the  equation.  The 
solution  is  expressed  as  a  flux  in  the  variables  z(cm), 
E(keV) ,  and  y  (direction  cosine).  The  front  surface 
solution  has  been  of  interest  in  this  program.  The  code 
has  considerable  input  in  the  form  of  atomic  parameters 
and  gridding  information,  the  latter  requiring  modifica¬ 
tions  from  run  to  run.  Based  on  the  difficulty  in 
selecting  good  electron  energy  grids,  one  would  have  to 
say  that  code  SXPJ>  is  not  a  user- oriented  code.  DNA  has 
recently  expressed  interest,  however,  in  having  the  code 
achieve  this  status  based  presumably  on  the  importance 
attached  to  soft  X-ray  photoemission,  the  vigorous  nature 
of  the  solution,  the  success  in  code  validation,  and  the 
many  material  models  developed  to  date.  A  significant 
step  was  made  during  the  past  year  toward  this  goal  by 
developing  an  algorithm  for  selecting  the  electron  energy 
grid.  This  proved  to  be  a  difficult  task  due  to  how 
highly  structured  the  photoelectron  source  spectrum  can 
be  and  constituted  a  major  effort  during  this  contract 


period.  A  description  of  the  algorithm  and  a  demonstra¬ 
tion  of  its  capability  will  be  given  in  this  report. 


The  above  discussions  cover  most  of  the  effort 
carried  out  under  the  given  contract.  There  were,  how¬ 
ever,  some  additional  investigations  reported  on  at  the 
1979  IEEE  meeting  on  nuclear  science.  One  of  these 
dealt  with  electron  backscatter  characteristics  for  At 
and  Au.  This  work  was  motivated  by  our  desire  to  fur¬ 
ther  critically  test  the  scattering  model  in  code  SXRP. 
Another  investigation  concerned  charge  and  energy  deposi¬ 
tion  in  Teflon  for  incident  soft  electron  sources.  Out 
of  this  work  has  come  a  material  model  for  teflon  and  a 
demonstration  of  code  SXRP's  capability  to  treat  external 
electron  as  well  as  X-ray  sources.  An  examination  was 
also  made  of  the  discrepancy  between  theory  and  experi¬ 
ment  for  PRS  radiation  incident  on  Ag.  Results  from 
these  assorted  studies  will  be  found  in  the  appendices 
which  contain  the  papers  as  published  in  the  IEEE  journal. 
From  recent  work,  we  have  also  generated  a  paper  in 
Journal  of  Applied  Physics.  A  copy  of  this  also  appears 
in  the  appendices. 
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Section  2 


FURTHER  CHARACTERIZATION  OF  SOFT  X-RAY 
PHOTOEMISSION  —  EMISSION  FOR  LINE  SOURCES 


In  this  section,  photoemission  properties  A£, 
Si(>2.  Au,  Ag  and  C  will  be  discussed  for  narrow 
X-ray  sources  effectively  simulating  line  sources. 
Selected  results  for  such  sources  have  previously  ap¬ 
peared  for  A£  and  Au^>^>^.  They  were  obtained  shortly 
after  the  respective  material  models  were  developed  to 
compare  with  available  published  data  in  the  form  of 
the  primary  yield  (electron/photon  for  electron  energies 
>  50  eV).  Here,  for  the  six  materials  listed  above,  we 
will  be  comparing  primary  yields  with  Bernstein's  mea¬ 
sured  values  considering  several  X-ray  lines  between 
1.26  and  5.4  keV.  Table  1  lists  the  lines  and  their 
energies. 


For  a  given  material,  results  have  not  been  ob¬ 
tained  for  all  lines.  We  had  planned  to  do  so  as  well 
as  at  other  energies  but  did  not  have  sufficient  com¬ 
puter  time  due  to  additional  investigations  related  to 
backscatter  from  A£  and  Au  and  charging  in  Teflon.  It 
is  our  desire  to  eventually  produce  a  soft  X-ray  photo¬ 
emission  handbook  which  will  contain  a  much  more  compre¬ 
hensive  set  of  results  of  the  type  to  follow. 

In  addition  to  primary  yields,  we  will  also 
be  presenting  photoemission  spectra  for  selected  lines. 
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Spectral  comparisons  with  experiment,  however,  will  not  be 
made  in  this  report.  Bernstein  has  been  measuring  photo¬ 
emission  spectra  and  some  preliminary  comparisons  have 
been  made.  The  reporting  of  such  information  shall  await 
his  own  presentation  of  the  measurements. 

It  can  be  stated,  however,  that  greater  varia¬ 
tions  appear  in  the  calculated  spectrum  which  cannot  be 
accounted  for  solely  by  differences  in  resolution.  In 
particular,  code  SXRP  gives  a  more  rapid  fall  off  with 
decreasing  electron  energy  away  from  spectral  maxima. 

This  difference  in  behavior  has  also  been  observed  with 
respect  to  a  limited  amount  of  data  by  Denisov  et  al.*^ 
(see  Strickland  et  al.®  for  details).  There,  as  here, 
we  observe  that  the  data  agree  more  closely  with  SXRP 
bulk  solutions.  The  cause  of  the  discrepancy  is  not 
understood  at  this  time. 

We  have  approximated  the  X-ray  line  source  by 
the  following  Gaussian  distribution: 

F(hv)  =  exp[-(hv  -  hv0)/a]2 

where  hvQ  is  the  line  energy  and  a  is  sufficiently  small 
that  the  finite  extent  of  F  has  little  effect  on  either 
the  yield  or  the  spectral  behavior  of  the  photoemission. 
Examples  of  source  spectra  will  be  given  in  the  next  sec¬ 
tion  which  illustrate  the  finite  extent  of  F. 

We  now  proceed  to  the  results  and  a  compari- 
son  with  the  measurements  of  Bernstein  and  Smithx  . 
Comparisons  with  other  data  may  be  found  in  two  previous 


« 


2  7 

reports  '  .  Figures  1-6  give  the  primary  yields  in  elec 
trons/photon  versus  photon  energy  for  At,  AA2O3,  SiC>2, 
Au,  Ag  and  C.  The  two  dashed  curves  come  from  a  simple 
empirical  model16*16  and  are  included  primarily  to  show 
where  absorption  edges  occur.  The  upper  of  the  two 
curves  implicitely  includes  electrons  down  to  10  eV 
while  the  lower  one  includes  them  down  to  only  100  eV. 
The  term  in  the  ertpirical  formula  giving  the  difference 
is  the  range.  Burke16  generated  his  results  using 
10  eV  ranges  and  obtained  favorable  agreement  with 
many  sets  of  primary  yield  data.  It  would  appear  to 
be  more  appropriate,  however,  to  use  the  100  eV  range 
since  100  eV  is  near  the  low  energy  cut  off  (by  conven¬ 
tion)  for  primary  electron  emission.  This  is  why  we 
include  the  empirical  results  for  the  100  eV  range  in 
the  figures.  We  note  that  the  two  curves  for  Au  are 
significantly  different  from  one  another  and  that  the 
more  appropriate  result  is  in  poor  agreement  with  the 
SXRP  results  and  measurements.  Based  on  such  observa¬ 
tions,  one  should  probably  not  take  too  seriously 
either  good  or  poor  agreement  between  the  empirical 
model  and  other  results,  at  least  at  the  lower  ener¬ 
gies  where  the  choice  of  cut  off  affecting  the  range 
becomes  important.  The  model  is  a  useful  tool,  however, 
for  estimating  the  general  behavior  of  the  primary 
yield  versus  photon  energy. 

Regarding  the  comparison  between  code  SXRP’s 
results  and  the  measurements,  we  view  the  extent  of 
agreement  to  be  satisfactory  considering  that  the  cal¬ 
culations  are  done  for  ideally  clean,  smooth  surfaces 
with  pure  material  beneath  them  which  can  only  approx¬ 
imate  experimental  conditions.  In  general,  our  yields 
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Figure  1  Primary  yield  for  A*  in  electrons/photon  versus 
photon  energy  In  keV.  BS  refers  to  the 
measurements  of  Bernstein  and  Smith.1  Dashed 
curves  come  from  empirical  models.  See  discus¬ 
sion  on  page  n. 
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Figure  2  Primary  yield  for  aluminum  oxide  In 
electrons/photon. 
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Figure  3  Primary  yield  for  silicon  dioxide  In 
electrons/photon. 
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Figure  6  Primary  yield  for  carbon  in  electrons/photon 
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fall  below  the  measurements  but  by  usually  less  than 
50%  (see  figures  for  Si02,  Au,  Ag  and  C).  Surprisingly, 
some  of  the  most  favorable  comparisons  have  been  made 
for  Au  which  is  the  most  difficult  to  model  of  those 
considered.  Experimentally,  however,  one  can  approach 
the  ideal  situation  (as  modelled)  more  closely  for  Au 
than  other  materials  such  as  At  which  rapidly  forms  an 
oxide  coating. 

The  largest  differences  are  found  for  Ag  be¬ 
low  its  L  edge  at  ^  3.5  keV.  Differences  as  large  as 
a  factor  of  two  exist  and  happen  to  occur  in  the  region 

* 

dominated  by  the  OWL  II'  PRS  X-ray  spectrum.  The 
OWL  II*  source  has  received  most  of  the  attention  both 
in  Bernstein's  and  our  program.  Previously,  direct  com¬ 
parisons  for  this  source  revealed  about  the  same  dis¬ 
crepancy  noted  above.  The  given  differences  have 
prompted  us  to  investigate  Ag  in  greater  detail.  Since 
good  agreement  had  been  previously  achieved  for  Au,  we 
decided  to  compare  transport  properties  between  Ag  and 
Au  by  allowing  for  the  same  electron  source  spectrum  in 
each  material.  We  found,  as  expected,  that  such  pro¬ 
perties  are  similar.  This  leads  one  to  the  conclusion 
that  differences  in  yield  for  these  materials  are  pri¬ 
marily  the  result  of  differences  in  respective  electron 
source  spectra.  Since  source  spectra  are  easy  to 
specify,  we  concluded  that  our  results  for  Ag  should 
be  as  accurate  as  those  for  Au.  Further  details  may 
be  found  in  Appendix  B  which  contains  figures  of  elec¬ 
tron  source  spectra  and  photoemission  spectra  used  in 
this  analysis. 
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We  now  turn  our  attention  to  spectral  charac¬ 
teristics  of  photoemission.  Results  will  follow  in  units 
of  electrons/photon-keV  versus  electron  energy  for  a  few 
select  line  sources.  The  purpose  of  presenting  such  in¬ 
formation  is  to  illustrate  the  complexity  which  can  exist 
in  a  photoemission  spectrum  and  the  diversity  possible  in 
going  from  one  X-ray  energy  to  another  and  from  one 
material  to  another.  Spectra  will  not  be  presented  for 
all  cases  treated  in  this  current  investigation.  They 
are  available  from  the  authors  and  will  probably  appear 
in  a  future  handbook  as  part  of  a  larger  set  of  results. 

The  selected  spectra  appear  in  Figures  7-13. 
Shown  in  each  figure  are  the  spectrum  and  its  corres¬ 
ponding  cumulative  yield  down  to  ^  .1  keV.  The  first 
two  spectra  shown  (Figures  7  and  8)  are  for  the  At 
K-line  (1.49  keV)  and  Cl  X-line  (2.64  keV)  incident  on 
At.  These  lines  straddle  the  At  K-edge  (1.56  keV)  and, 
as  can  be  seen,  produce  significantly  different  spectra 
and  yields.  The  features  are  identified  in  each  figure 
and  are  either  photoelectron  (single  letter  for  shell 
from  which  electron  is  ejected)  or  Auger  (three  letters). 
The  first  letter  in  the  Auger  designation  gives  the 
shell  from  which  a  photoelectron  is  ejected.  The  second 
letter  gives  the  shell  from  which  an  electron  originates 
that  drops  into  the  created  vacancy.  The  third  letter 
gives  the  shell  from  which  the  Auger  electron  is 
ejected. 

A  demonstration  of  material  effect  on  the 
photoemission  spectrum  is  provided  by  Figure  8  for  At 
just  discussed  and  Figures  9  and  10  for  A&203  and  Si02, 
all  of  which  refer  to  the  Cl  K-line  source.  Fairly 
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gure  9  Photoemission  spectrum  and  its  cumulative  yield 
for  a  2.64  keV  Gaussian  photon  source  Incident 
on  AtgOj* 
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Figure  11  Photoemission  spectrum  and  its  cumulative 

yield  for  a  4.12  keV  Gaussian  photon  source 
incident  on  Au. 
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Figure  13  Photoemission  spectrum  and  Its  cumulative 

yield  for  a  1.49  keV  Gaussian  photon  source 
Incident  on  C. 


complex  spectra  are  produced  from  the  oxides  due  to  the 
added  features  from  oxygen.  Overall,  the  presence  of 
oxygen  reduces  the  strength  of  photoemission  since  oxy¬ 
gen  does  not  tend  to  photoionize  in  and  above  the  soft 
X-ray  region  as  readily  as  higher  Z  materials  such  as 
At. 

Figures  11  and  12  show  spectra  for  the  Sc 
K-line  incident  on  Au  and  Ag.  Again,  such  results 
allow  for  a  demonstration  of  material  effect,  this 
time  for  much  high  Z  materials  (Z^u  *  79;  Z^g  *  47). 
Finally,  Figure  13  provides  an  example  of  photoemission 
from  C.  The  chosen  source  is  the  A £  K-line  at  1.49  keV. 
Unlike  the  higher  Z  materials,  which  can  produce  sig¬ 
nificantly  different  emission  spectra  as  hv  varies  over 
the  soft  X-ray  region,  C  will  only  produce  the  three 
features  shown  with  K-photoelectrons  dominating  the 
emission. 


Section  3 


INCORPORATING  A  SELF-CONTAINED  ENERGY  GRID 
SELECTION  SCHEME  INTO  CODE  SXRP 


One  of  the  more  difficult  tasks  undertaken 
in  our  overall  program  will  be  discussed  in  this  sec¬ 
tion.  That  task  has  been  to  develop  a  scheme  which 
provides  a  satisfactory  electron  energy  grid  for  trans¬ 
port  calculations  using  code  SXRP.  We  can  report  that 
we  have  successfully  developed  such  a  scheme.  As  a  re¬ 
sult  ,  considerable  time  has  been  and  will  be  saved  in 
preparing  the  input  data  for  SXRP  runs.  The  scheme 
provides  good  energy  grids  for  either  narrow  or  broad 
photon  and  electron  sources  for  all  materials  so  far 
modeled. 


The  difficulties  in  developing  the  above 
scheme  stem  from  the  following  conditions: 

•  limited  number  of  available  grid 
points  (currently  40) 

•  complex  structure  of  photoelectron 
source  spectra,  and 

•  requirement  that  changes  in  grid 
interval  sizes  be  gradual. 

Figures  14  and  15  show  examples  of  the  complexity  one 
can  encounter  in  the  source  spectrum.  For  the  case 
considered  in  Figure  14,  a  5.4  keV  Gaussian  X-ray 
source  is  incident  on  Au.  Figure  15  shows  the  spectrum 
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for  a  1.49  keV  Gaussian  source  incident  on  Ag.  The 
grid  points  shown  were  selected  by  the  code.  The  points 
on  axis  at  the  bottom  of  the  figures  are  actually  zeros. 
For  completeness,  the  corresponding  photoemission  spectra 
are  shown  in  Figures  16  and  17. 

The  third  condition  above  is  a  requirement 
general  to  most  transport  formulations.  Specific  to 
that  of  SXRP,  it  centers  about  the  representations  of 
the  solution  <J>  (z,E,y)  in  the  Boltzmann  collision  inte¬ 
gral.  To  cast  the  integral  into  matrix  form,  <p  is 
given  quadratic  dependence  versus  tnE  and  linear  depen¬ 
dence  versus  direction  cosine  y  over  each  E,y  cell. 

With  this  quadratic  dependence,  one  must  be  careful  not 
to  abruptly  change  the  grid  spacing  in  energy.  Other¬ 
wise,  unphysical  structure  may  appear  in  the  solution. 

We  will  now  proceed  to  describe  the  grid 
selection  scheme  in  two  stages.  A  brief  overview  will 
first  be  given  in  this  paragraph  then  to  be  followed 
by  selected  details.  We  say  selected  since  an  attempt 
will  not  be  made  to  provide  a  detailed  documentation. 

Such  an  effort  would  be  difficult  due  to  the  degree  of 
logic  required  in  the  scheme.  The  high  level  called 
for  arises  from  internal  decisions  such  as  when  to 
merge  features  due  to  close  proximity  and  from  fulfill¬ 
ing  such  conditions  as  gradual  grid  interval  variation 
throughout  the  entire  energy  range.  We  begin  by  noting 
that  the  scheme  works  on  a  logarithmic  scale  in  energy. 
This  is  due  to  the  way  the  solution  $  is  represented  in 
the  collision  Integral  (as  noted  above)  and  the  capa¬ 
bility  of  code  SXRP  to  handle  energy  ranges  spanning 
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several  decades.  Grid  selection  on  a  logarithmic  scale 
allows  naturally  for  an  increase  in  resolutions  with 
decreasing  energy.  The  grid  selection  begins  with  spe¬ 
cification  of  what  shall  be  referred  to  as  principal 
points.  These  points  give  the  energies  at  which  the 
Auger  and  photoelectron  features  peak.  The  initial  set 
of  points  is  then  examined  in  terms  of  strength  and 
proximity  of  the  various  features.  If  the  strength  of 
a  given  feature  is  less  than  a  given  level  relative  to 
the  other  features,  its  principal  point  is  removed. 

If  two  or  more  features  are  determined  to  be  close  to¬ 
gether  based  on  both  the  amount  of  overlap  and  their 
separation  relative  to  the  full  energy  range,  all  ex¬ 
cept  one  principal  point  are  removed.  A  set  of  grid 
points  is  then  specified  for  each  principal  point  to 
characterize  the  given  feature.  Five  such  points, 
equally  spaced,  are  used  for  Auger  features.  The  num¬ 
ber  of  points  for  photo  features  is  variable  but  is  no 
less  than  three.  If  a  given  photo  feature  is  sharply 
peaked,  more  points  are  inserted  within  the  structure. 
The  final  step  is  to  add  points  between  the  adjacent 
principal  point  subsets  in  a  way  that  achieves  a  smooth 
variation  in  adjacent  grid  spacings.  We  currently  do 
not  allow  the  ratio  of  adjacent  intervals  to  lie  out¬ 
side  the  range  .5  to  2.  To  achieve  the  required 
variations  and  keep  the  total  number  of  grid  points 
less  than  or  equal  to  40,  an  iteration  technique  is 
introduced  which  allows  for  adjustments  in  both  the 
principal  point  subsets  and  the  points  between  them. 


The  mathematical  details  of  the  essential 
logic  for  grid  selection  will  be  given  in  this 
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paragraph.  Since  the  adjacent  principal  points  are 
well-separated,  we  can  assume  that  the  interval  be¬ 
tween  them  is  large  enough  that  principal  point  sub¬ 
sets  do  not  overlap.  We  want  to  add  points  starting 
from  the  sides  of  the  interval.  A  comparison  is 
first  made  between  the  sizes  of  the  existing  spacings 
on  both  sides.  A  point  is  added  at  the  side  adjacent 
to  the  smaller  spacing  in  such  a  way  that  the  size  of 
the  new  spacing  is  twice  as  large  as  the  previous  one 
on  this  side.  A  comparison  is  again  made  between  the 
size  of  this  new  spacing  and  the  one  on  the  other  side, 
with  a  point  then  added  to  the  side  of  the  smaller 
spacing  in  the  same  fashion  as  before.  This  insertion 
procedure  is  stopped  when  the  latest  point  from  one 
side  overshoots  a  point  on  the  other  side.  For  ex¬ 
ample,  in  the  diagram  below,  point  d  from  the  right 
hand  side 


overshoots  the  point  D  inserted  from  the  left  hand  side. 
The  logic  requires  the  removal  of  the  two  points  D  and 
d  and  determination  of  a  new  grid  point  at  D'  such  that 


X  =  CD' 


AXj+AXg 


as  shown  at  the  top  of  the  next  page. 
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*sAx  j 

i- 

i-h' 

I- - 1 

The  three  spacings  ratios  CD":BC,  D"C:CD", 
and  cb:D"C  are  between  .5  and  2.  To  prove  this,  we 
make  use  of  the  two  given  inequalities 

Ax1  <  Ax2  <  2Ax^ 


and 


2Ax2  <  SL  <  AXj^  +  Ax2 


Then 


CD"  _  x  =  2d  <  2 

BC  4Ax1  ~  Ax^+AXg 


On  the  other  hand 

C"D"  =  2d  4AX1  =  4  >  4 

BC  AXj^+AXg  Axj+AXg  Axg  3 

1+Ax^ 


For  the  next  ratio, 


rc  =  d^x  =  Ax2 

CD"  x  ~  Ax^ 
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which  is  between  1  and  2.  The  last  ratio  cb/D'c  is, 
due  to  the  way  x  is  defined,  equal  to  the  first  ratio 
CD'/BC. 
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Section  4 


ADDITIONAL  STUDIES  DOCUMENTED  IN  PUBLICATIONS 


The  purpose  of  this  section  is  to  introduce 
and  summarize  the  three  papers  represented  in 

P 

Appendices  A-C.  The  first  paper  by  Strickland  and  Lin 
demonstrates  the  capability  of  code  SXRP  to  treat  ex¬ 
ternal  soft  electron  sources.  The  two  topics  addressed 
are  primary  backscatter  from  Ai  and  Au  and  charge  and 
energy  deposition  in  Teflon.  The  motivation  for  the 
backscatter  analysis  comes  from  our  desire  to  critically 
test  our  scattering  model  for  materials  of  high  and  low 
atomic  number.  The  incident  energy  range  considered  is 
from  .5  to  10  keV  and  good  agreement  is  achieved  in  com¬ 
parison  with  data.  One  comparison  in  the  backscattered 
spectrum  is  made  with  code  POEM1-7  and  there,  good  agree¬ 
ment  is  also  obtained.  The  motivation  for  the  calcula¬ 
tions  of  deposition  in  Teflon  comes  from  a  charging 
analysis  by  Beers  et  al.*®.  A  simple  material  model  for 
Teflon  is  first  developed  but  is  not  reported  on  in  the 
paper.  This  will  await  further  refinement  which  we 
hope  to  do  as  part  of  a  future  investigation  of  electron 
transport  in  polymers.  Charge  and  energy  deposition 
profiles  are  generated  with  code  SXRP  for  2.5  and  5.0  keV 
electrons  incident  on  Teflon  and  compared  with  the  POEM 
results  used  by  Beers  et  al.  The  main  interest  is  in 
differences  in  penetration  characteristics  since  these 
have  a  bearing  on  charge  buildup.  As  expected,  the 
SXRP  profiles  extend  further  into  the  material  due  to 


the  code's  treatment  of  straggling  but  the  exhibited 
differences  are  not  large. 


The  second  paper  by  Lin  and  Strickland0  dis¬ 
cusses  photoemission  from  Ag  and  C  for  the  OWL  II " 
source.  Spectral  and  primary  yield  comparisons  are 
made  with  the  measurements  by  Bernstein*!.  These 
comparisons  were  previously  reported  in  a  recent  DNA 
final  report7.  The  extended  analysis  of  Ag  is  re¬ 
ported  for  the  first  time  in  the  parer.  It  came  about 
from  the  factor  of  two  difference  in  primary  yield  be¬ 
tween  the  calculation  and  measurement.  Following  the 
analysis,  we  better  understand  transport  properties  in 
Ag  and  have  more  confidence  in  the  calculations.  The 
difference  noted  above  still  remains  to  be  determined. 

The  third  paper  by  Lin  and  Strickland10  re¬ 
ports  our  latest  material  model  for  Ag.  The  discussion 
here  is  more  theoretical  than  in  the  report  by 
Strickland  and  Lin  which  introduced  the  model7.  The 
model  parameters  in  the  paper  also  extend  into  the 
very  soft  X-ray  region. 
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INTRODUCTION 

In  this  paper,  we  will  address  the  basic  problea 
of  low  energy  electron  transport  in  materials.  By  low 
energy,  we  are  here  referring  to  the  low  to  subkilo¬ 
volt  energy  region. 

As  a  source,  we  consider  electrons  Incident  with 
distributions  sufficiently  narrow  in  energy  to  ade¬ 
quately  simulate  a  mono-energetic  condition.  For  the 
external  electron  source,  three  basic  transport  quan¬ 
tities  of  general  Interest  are 

#  the  primary  backscatter  yield  (energy  dif¬ 
ferential  and  energy  integrated) 

e  primary  current  profile,  and 

e  dose  or  energy  deposition  profile. 

For  selected  materials  we  will  be  addressing  each  of 
these  quantities.  To  do  so,  we  apply  the  transport 
code  SXRP  which  provides  for  a  rigorous  description 
of  transport  down  to  sub-keV  energies.  From  SXRP 
will  come  backscatter  results  for  Ai  and  Au  and  cur- 
;nt  and  dose  profile  information  for  Teflon  (CF^). 

The  code  SXRP  and  its  applications  are  by  now 
fairly  well  documented, Its  primary  application 
to  date  has  been  to  the  problem  of  soft  X-ray  photo- 
emission.  The  transport  description  is  provided  by 
the  Boltzmann  equation  which  is  solved  for  the  elec¬ 
tron  flux  <p(z,E,D)  in  e/cm^-s-eV-sv.  The  arguments 
of  0  are  depth  (z),  energy  (E)  and  direction  cosine 
(U) .’  The  given  formulation  allows  for  the  produc¬ 
tion  of  photo-.  Auger  and  Compton  electrons  by 
photons  and  for  the  following  effects  involving 
electrons: 

a  scattering  through  any  angle, 

s  discrete  energy  loss,  and 

s  production  of  Auger  and  secondary  electrons. 

SXRP  is  thus  capable  of  investigating  transport 
effects  down  to  energies  well  below  the  cutoff 
imposed  by  the  continuous  slowing  down  approximation 
(CSDA)  which  is  commonly  used  in  other  formulations. 

There  are  two  prime  motivating  factors  for  the 
work  to  be  discussed  in  this  paper: 

(1)  electron  backscatter  is  sensitive  to  the 
scattering  and  energy  loss  characteristics 
of  a  given  material  and  consequently  pro¬ 
vides  a  good  test  on  the  material  model 
needed  to  perform  the  transport,  and 


The  backscatter  analyaia  la  thus  basically  a  material 
model  validation  exercise.  AC  and  Au  have  been  chosen 
for  this  analysis  since  most  of  our  efforts  to  date 
have  concentrated  on  these  materials.  Beyond  validat¬ 
ing  our  models,  the  analysis  also  improves  our  basic 
understanding  of  transport  properties  down  to  sub-keV 
energies  and  demonstrates  the  applicability  of  code 
SXRP  to  external  electron  at  well  as  X-ray  sources. 

The  work  by  Beers,  et  al.  addressed  the  general 
problea  of  spacecraft  charging.  One  part  of  their 
analysis  dealt  with  the  equilibrium  charge  distribu¬ 
tion  in  Teflon  for  a  2.5  keV  electron  source.  The 
needed  current  and  dose  profiles  were  obtained  using 
the  Monte  Carlo  transport  code  POEM.’  Since  the  code 
SXRP  is  well  suited  to  treating  a  source  as  low  in 
energy  as  the  above,  a  simple  Teflon  model  was  devel¬ 
oped  and  applied  to  the  same  problea.  Me  were  pri¬ 
marily  Interested  in  seeing  how  different  the  penetra¬ 
tion  characteristics  would  be  between  the  tiro  approaches. 
This  work  will  report  our  first  such  results  and  show  a 
comparison  with  the  POEM  results. 


MATERIAL  MODEL  PARAMETERS 

The  following  parameters  are  needed  in  our  formu¬ 
lation  to  describe  electron  interactions: 

e  elastic  scattering  differential  lnverae 
mean  free  path  (DIMFP), 

e  plasmon  DIMFP, 

e  conduct ion/ valance  band  DIMFP,  and 

e  inner  shell  DIMFP's. 

A  discussion  will  follow  which  describes  how  we 
currently  specify  the  elastic  DIMFP.  We  have  previ¬ 
ously  described  how  the  other  parameters  are  specified 
and  thus  their  description  here  will  be  brief.  Basic¬ 
ally,  the  plasmon  and  conductlon/valance  band  DIMFP's 
are  obtained  from  the  dielectric  response  function  as 
modeled  by  Ritchie  and  his  colleagues.  11  The  inner 
shell  IMFP's  (DIMFP's  integrated  over  either  energy 
loss  or  secondary  electron  energy)  are  obtained  from 
their  corresponding  photolonixatlon  cross  sections  by 
the  method  described  by  Lin,  et  ai.*  The  differential 
in  energy  loss  is  that  given  by  Strickland, 
The  results  for  At  and  Au  from  applying  the 
techniques  have  already  been  presented  elsewhere.  •’* 
The  Teflon  IMFP's,  however,  arc  new  but  relatively 
untested  and  we  thus  defer  to  a  later  date  their  pre¬ 
sentation. 

We  use  the  screened  Rutherford  formula  to  s,  ify 
the  elastic  scattering  DIMFP.  Its  form  is  given  v 


(2)  the  previous  work  by  Beers,  et  al.**  on 
electron  charging  in  Teflon. 


Xeias^®'®)  ‘  *elas(E)  p(E'0) 


(D 
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where 


p(E,9>  -  - i - -  (2) 

(i-co»0+2n) 

p(E,6)  ia  the  normalized  Rutherford  formula  and  n  la 
the  acreenlng  parameter  which  contalna  the  energy  de¬ 
pendence  In  p.  Keia>(E)  haa  the  form 


,2  4 

Z  a  a 

v2p2  n<nfi) 


(3) 


where  v  and  p  are  reapectlvely  the  electron  velocity 
and  momentum  corrcapondlng  to  E,  Z  ia  the  atonic  num¬ 
ber,  and  n  ia  the  material  number  density.  The 
screening  parameter  n  may  be  represented  in  the  form 


r>(e) 


4.3 


,2/3  % 
Z  E 


(4) 


where  E  ia  in  eV  and  nc  is  a  function  of  the  screening 
potential.  The  screening  potential  in  conducting 
solids  can  be  adequately  represented  by  a  single 
Yukawa  potential  which  leads  to  constant  n..*  *  We 
have  previously  used  nc  -  1  since  values  close  to 
unity  have  been  applied  before. 14,15  We  now  use  the 
value  3.2,  however,  since  it  is  based  directly  on  con¬ 
siderations  of  the  screening  potential  itself. * 3 

The  elastic  IMFP  from  Eq.  3  and  the  screening 
parameter  n  are  shown  in  Fig.  1  for  both  At  and  Au 
over  the  range  of  0.1  to  10  keV.  The  point  of  showing 
these  results  Is  to  emphasize  the  contrast  between  At 
and  Au  in  their  scattering  characteristics.  Both  the 
elastic  collision  frequency  and  mean  angle  of  scatter¬ 
ing  per  collision  are  larger  for  Au.  We  have  not 
bothered  showing  Keiag  and  n  for  Teflon  since  their 
variations  with  energy  will  look  similar  to  those  in 
Fig.  1.  The  needed  results  are  obtained  with  the 
above  equations  for  Z  values  of  6(C)  and  9(F)  and  for 
the  appropriate  number  densities  of  these  elements. 

A  word  of  caution  is  appropriate  here  regarding 
these  results  at  low  energies.  Eq.  (3)  with  n  given 
by  Eq.  (4)  may  become  increasingly  inaccurate  with 
decreasing  energy  since  these  expressions  apply  to  the 
Born  region  which,  though  Z  dependent,  is  generally 
above  the  low  keV  region.  For  lack  of  definitive  and 
comprehensive  information  on  Keias  at  low  energy,  we 
choose  to  stay  with  expressions  (2)-(4)  which  not  only 
are  particularly  convenient  to  work  with,  but  also  ex¬ 
press  the  isotropic  nature  of  acoustic  phonon  scatter¬ 
ing  generally  accepted  for  electron  energies  of  a  few 
eV. 


maximum  number  of  angles  which  SXRP  can  handle,  we 
choose,  for  the  sake  of  saving  computing  time,  instead 
to  apply  corrections  to  our  backscatter  results  for 
comparisons  with  measurements.  The  corrections  are 
based  on  available  data  and  will  be  further  discussed 
below. 

As  distributions  finite  in  width,  we  have  chosen 
normalized  Gaussian  distributions  for  characterizing 
the  incident  electron  flux.  Designated  by  the 
incident  flux  is  given  by 

♦0(E,h)  •  f (E)g(p)  el/cm2-s-eV-sr  (5) 

where 

"E>  -  TO  *  “  <•> 

and 

(1-u)2 

8<u)  "  6 (tt-6)  e  (7) 

These  functions  are  normalized  to  give  a  total  incident 
flux  of  1  el/cm^-s.  The  parameters  a  and  6  control  the 
widths  of  the  distributions  while  the  parameter  E0  gives 
the  energy  at  which  peaks.  We  have  considered  a  num¬ 
ber  of  sets  of  E0  and  a  which  are  shown  in  Table  1 .  The 
absolute  width  of  f (E)  is  allowed  to  broaden  as  E0  in¬ 
creases  which  leads  to  the  variation  in  a  as  shown. 

We  consider  only  one  3  value  In  this  work,  namely, 
B*0.2.  We  feel  that  this  is  about  as  small  as  the 
value  should  be,  given  a  maximum  of  ten  angles  per 
hemisphere.  To  adjust  our  yields  to  normal  incidence, 
we  fold  the  following  expression  into  the  above  g(u): 

Y(U)  -  B  (8) 

This  expression  comes  from  Darlington1*  who  assigned  B 
the  value  0.891.  Darlington  arrived  at  Expression  (8) 
and  the  above  value  of  B  by  considering  data  for  inci¬ 
dent  energies  above  9.3  keV.  We  expect  them  to  remain 
adequate  down  to  at  least  1  VeV  since  yield  data,  in 
general,  do  not  show  much  variation  from  the  10  keV 
range  on  down  to  this  energy.  Departures  that  may  take 
place  should  be  In  the  direction  of  weaker  p  dependence 
since  the  mean  scattering  angle  Increases  with  decreas¬ 
ing  energy. 


ENERGY  AND  ANGULAR  DEPENDENCE 
OF  THE  INCIDENT  ELECTRON  FLUX 

We  wish  to  consider  incident  electron  fluxes  nar¬ 
row  in  both  energy  and  angle  since  in  the  experimental 
situation,  nearly  monoenergetic  and  normal  incidence 
conditions  usually  exist.  SXRP,  however,  can  only  ap¬ 
proximate  these  conditions,  the  extent  of  which  depends 
on  the  number  of  grid  points  chosen  in  the  above  varia¬ 
bles.  The  code  is  currently  set  up  to  handle  a  maximum 
of  40  energy  and  20  angular  points.  For  past  photo- 
emission  work,  these  numbers  have  been  more  than  ade¬ 
quate  for  satisfactory  comparison  of  results  with  mea¬ 
surements.  The  same  still  holds  true  for  backscatter 
results  in  terms  of  the  energy  variable  since  back¬ 
scatter  characteristics  are  slowly  varying  in  energy. 
This  is  not  the  case,  however,  in  the  angular  variable 
due  to  the  sensitivity  of  the  backscatter  yield  with 
the  angle  of  incidence.  Rather  than  Increase  the 


BACKSCATTER  FROH  At  AND  Au 

The  information  to  be  presented  in  this  section, 
although  not  contributing  significant  new  insights  into 
emission  characteristics  of  materials  for  external  elec¬ 
tron  sources,  is  important  for  assessing  the  overall 
validity  of  our  formulation,  material  models,  and  the 
SXRP  code  itself.  A  comparison  of  backscatter  results 
with  measurements  provides  one  of  the  most  critical 
tests  that  can  be  made  on  the  scattering  model  in  a 
given  formulation. 

The  significance  of  such  a  test  lies  in  the  fact 
that  scattering  is  one  of  the  most  Important  processes 
in  electron  transport  be  it  either  to  problems  dealing 
with  photon  or  electron  sources.  One  cannot  begin  to 
reasonably  predict  current  and  dose  characteristics 
within  materials  for  electron  sources  without  taking 
it  into  consideration.  The  situation  is  less  critical 
for  photon  sources  but  scattering  still  does  have  a 
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pronounced  affect  on  such  quantities  as  prtaary  yield 
end  lta  aaaoclated  spectral  behavior. 5 

Hera,  vs  will  concentrate  on  primary  yield  results 
for  At  and  Au  between  0.5  and  20  keV.  This  range  la 
and  has  been  of  particular  interest  to  ua  in  character¬ 
ising  photoaadssion  properties  of  aaterlala  such  as  the 
above  aa  veil  as  AI2O3,  S3O2.  Ag,  Cu  and  C.  For  code 
validation  purposes,  At  and  Au  are  particularly  desir¬ 
able  since  they  provide  cases  of  low  Z  (ZAj  •  11)  and 
high  Z  (Zgy  ■  79)  and  are  aaterlala  for  which  exten¬ 
sive  aeaaureaenta  have  been  nade. 

Fig.  2  presents  our  calculated  primary  backscatter 
yield  curves  for  At  and  Au  fron  O.S  to  20  keV.  They 
appear  together  with  the  aeasureaents  taken  fron 
Fitting17  and  with  predictions  fron  the  code  POEM. 

Other  aeasureaents  have  also  been  reported  in  the  lit¬ 
erature  but  are  not  included  here  since  nost  of  then 
closely  follow  Fitting's  results  (see  Burke1®  for  s 
recent  tabulation  of  the  data) .  The  calculated  values 
include  contributions  down  to  0.1  keV  and  have  been  ad¬ 
justed  to  normal  Incidence  using  Eq.  (8)  with  Y(l)  in 
that  equation  taken  fron  the  data  in  Fig.  2.  POEM 
yield  results  for  Au  are  not  shown  below  10  keV  and 
for  Ai  not  below  2  keV  since  they  begin  to  rapidly 
decrease  below  these  energies.  The  lower  Unit  of 
applicability  of  this  code  is  seen  to  have  a  pronounced 
Z  dependence.  This  indicated  trend  would  suggest,  how¬ 
ever,  that  POEM  should  give  an  adequate  transport 
description  down  to  a  few  keV  for  the  third  material 
addressed  in  this  paper,  naaely.  Teflon.  We  will  show 
that  this  is,  in  fact,  the  case  in  the  next  section. 
Returning  to  Fig.  2,  the  overall  agreeaent  between  the 
SXRF  yields  and  the  data  is  good  although  trends  with 
energy  are  soaevhat  different.  If  one  accepts  the  data 
trends,  the  coaparlsons  suggest  we  overestimate  scatter¬ 
ing  for  Au  and  underestimate  it  for  Ai  at  low  energies 
relative  to  the  loss  processes.  If  this  is,  in  fact, 
the  case,  we  do  not  know  at  this  time  whether  our 
scattering  description  or  our  energy  loss  description 
is  the  more  responsible  for  the  given  differences. 

Figure  3  gives  examples  of  differential  yields  or 
backacattered  spectra  for  5  keV  electrons  incident  on 
the  two  aaterlala.  Cuaailatlve  yields  are  shown  below 
the  spectra  and  give  the  unadjusted  primary  yield 
values  on  the  left.  The  Au  apectrua  is  seen  to  be 
less  energy  degraded  due  to  the  stronger  scattering 
properties  of  this  material. 

Figure  A  offers  a  comparison  between  SXRP  and  POEM 
spectra  for  the  sane  5  keV  source  incident  on  Ai.  Both 
sets  of  results  apply  to  the  non-normal  Incidence  case 
described  by  g(p)  in  Eq.  (7).  Good  agreeaent  in  spec¬ 
tra  and  yields  is  obtained  in  spite  of  the  diverse 
methods  applied  and  strongly  energy  degraded  nature 
of  the  spectra. 


CURRENT  AMP  POSE  PROFILES  IN  TEFLON 

In  this  section,  we  consider  2.3  and  3  keV  elec¬ 
trons  incident  on  Teflon  which  wo  have  simply  modeled 
as  a  alxture  of  C  and  F  in  their  appropriate  amounts. 
Calculated  current  and  dose  profiles  will  be  presented 
and  coaparad  with  similar  results  obtained  with  the 
POEM  code.  Our  initial  interest  in  using  SXRP  to  ad¬ 
dress  this  problem  was  to  see  how  much  affect  straggling 
would  have  on  these  profiles.  POEM  does  not  treat  this 
affect  and,  consequently,  its  calculated  currant  and 
dose  become  identically  aero  at  a  depth  of  one  electron 
range  into  tha  material.  SXRP  la  ideally  suited  to  In¬ 
vestigate  depth  effects  since  it  rigorously  allows  for 
straggling  through  its  discrete  energy  loss  description 
and  because  of  its  method  of  solving  the  Bclttaarm 
transport  aquation.  The  method  is  based  on  a 


eigenvalue  approach  Which  can  be  expected  to  be  super¬ 
ior  to  either  finite  difference  or  Monte  Carlo  schemes 
for  specifying  transport  characteristics  deep  in  the 
material. 

The  rasulta  to  follow  are  based  on  two  simple  con¬ 
straints  we  have  forced  upon  our  initial  material  modal 
for  Taflon.  These  are: 

s  the  inelastic  DIMFP's  give  the  Bethe  stopping 
power  above  several  keV,  and 

a  the  elastic  DIMFP,  for  the  given  inelastic 
DIMFP's,  give  a  good  primary  backscatter 
yield  (A.13Z  in  low  keV  region). 

These  are  iaportant  constraints  since  by  experience,  we 
have  found  that  penetration  characteristics  are  quite 
sensitive  to  the  ratio  of,  say,  the  elastic  IMFP  to  the 
stopping  power. 

Figure  3  shows  SXRP  and  POEM  current  and  dose  pro¬ 
files  for  2.3  keV  electrons  incident  on  Teflon.  Figure 
6  shows  similar  results  for  5  keV  electrons.  In  general, 
good  agreeaent  is  obtained  between  the  two  different 
transport  models.  The  SXRP  results  are  seen  to  extend 
deeper  into  the  aaterlal,  aa  expected,  due  to  the  effect 
of  straggling.  A  pronounced  difference  is  not  produced, 
however,  by  this  effect.  The  largest  differences  occur 
in  dose  near  the  front  face.  Some  enhancement  in  the 
SXRP  dose  will  occur  due  to  its  source  not  being  purely 
normally  incident.  This  will  not,  however,  explain  the 
bulk  of  the  difference.  Ue  believe  it  is  due  to  energy 
deposition  below  M  keV  which  is  Included  in  SXRP's  de¬ 
scription,  but  not  in  POEM's  due  to  the  high  energy  na¬ 
ture  of  that  code. 

To  suaaarlze  this  section,  overall  good  agreement 
between  SXRP  and  POEM  has  been  obtained  in  current  and 
dose  profiler  in  Teflon.  We  wish  to  stress,  however, 
that  our  material  model  for  Teflon  has  just  been 
developed,  is  staple  in  nature,  and  as  of  now  is  not 
well  validated.  Furthermore,  we  have  found  that  the 
current  and  dose  profile  behaviors  are  quite  sensitive 
to  changes  in  the  scattering  and  energy  loss  descrip¬ 
tions  within  this  model.  Thus,  as  our  model  becomes 
more  refined  and,  in  turn,  better  validated,  some 
changes  may  occur  in  our  results. 


SUMMARY  AND  CONCLUSIONS 

In  this  paper,  we  have  reported  on  two  recent 
analyses  carried  out  using  the  code  SXRP  to  treat 
external  electron  sources.  The  first  analysis  con¬ 
cerned  backscatter  from  Ai  and  Au  and  was  undertaken 
to  critically  test  the  scattering  and  energy  loss 
characteristics  in  our  aaterlal  models.  We  feel  the 
test  was  successful  based  on  the  agreement  obtained 
with  published  primary  backscatter  yield  data  between 
0.3  and  20  keV.  Since  photoealssion  characteristics 
are  less  sensitive  to  transport  parameters  than  those 
for  backscatter,  these  tested  models,  as  used  by  code 
SXRP,  should  be  especially  valid  for  the  problem  of 
photoealssion .  This  has  already  been  indicated  by  the 
agreeaent  obtained  with  the  available  published  photo- 
emission  data. 

The  second  analysis  concerned  current  end  dose  pro¬ 
files  in  Teflon  and  was  undertaken  to  test  one  aspect  of 
a  charging  calculation  by  Beers,  et  al.®  concerned  with 
the  larger  problem  of  spacecraft  charging.  The  points 
of  this  snalysls  were  to  determine  the  applicability  of 
code  SXRP  to  specifying  Internal  transport  characteris¬ 
tics,  but  more  specifically  to  see  how  iaportant 
straggling  is  to  the  above  profiles  since  their  fall- 
off  behavior  deep  in  the  material  is  iaportant  to  the 
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applied  charging  aodal.  Baara,  at  al.  had  uaad  the 
t rasa port  coda  POM  to  obtain  their  currant  and  doee 
prof Ilea  which  doaa  not  permit  an  electron  to  go  beyond 
one  CSOA  range.  Ha  found  that  the  SXRP  calculated  p ro¬ 
ll  lea  do  extend  deeper  into  the  material  due  to  strag¬ 
gling,  but  that,  in  general,  the  agreement  with  the 
FOB!  raaulta  was  good  except  near  the  front  surface. 
Thia  conclusion  is  somewhat  conditional  basad  on  our 
present  material  model  for  Teflon  which  has  just  been 
developed  and  is  in  need  of  further  work. 
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TABLE  1.  ENERGY  PARAMETERS  FOR 

INCIDENT  ELECTRON  FLUXES 


0.1  1  10 


ENERGY  <kcv> 

FIGURE  1.  ELASTIC  INVERSE  MEAN  FREE  PATHS  AND 
SCREENING  PARAMETER  n  FOR  At  AND  Au 
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SCREENING  PARAMETER 


ENERGY  <koV>  ENERGY  tkeV) 


FIGURE  2.  ELECTRON  BACKSCATTER  YIELDS  FOR 
A£  A NO  Au  FROM  0.5  TO  20  keV 


ENERGY  (koV) 


FIGURE  3.  DIFFERENTIAL  AND  CUMULATIVE  BACKSCATTER 
YIELDS  AS  CALCULATED  BY  CODE  SXRP  FOR  A 
5  keV  GAUSSIAN  ELECTRON  SOURCE  INCIDENT 
ON  A4  AND  Au 

These  results  have  not  been  adjusted  to  normal 
Incidence  from  the  Gaussian  type  of  incident 
angular  dependence. 


FIGURE  4.  DIFFERENTIAL  YIELDS  AS  CALCULATED  BY  THE 
CODES  SXRP  AND  POEFI  FOR  A  5  keV  GAUSSIAN 
ELECTRON  SOURCE  INCIDENT  ON  A£ 

Also  shown  are  SXRP's  cumulative  yield  and  POEM's 
primary  yield.  These  results  (including  POEM's), 
like  those  in  Figure  3,  have  not  been  adjusted  to 
normal  Incidence. 


REPTH  (urn) 

FIGURE  5.  CURRENT  AND  DOSE  PROFILES  AS 
CALCULATED  BY  THE  CODES  SXRP 
AND  POEM  FOR  2.5  keV  ELECTRONS 
INCIDENT  ON  TEFLON 

The  dots  give  the  POEM  results  while  the  dashed 
curves  give  fits  to  these  results. 
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PHOTOEMISSION  FROM  Ac  AND  C  FOR  AN 
EXPLODING  NIRE  RADIATOR  SOURCE* 


D.  L.  Lin  nnd  D.  J.  Strickland 
RADIATION  AND  ELECTR0NA6NETICS  DIVISION 
SCIENCE  APPLICATIONS.  INC. 

V ienna ,  Virginia  22180 


INTRODUCTION 

In  chi*  paper,  we  will  discus*  the  problem  of 
sofc  X-ray  photoemlsslon  for  Che  material*  Ag  and  C. 

In  particular,  we  will  be  presenting  calculated  photo- 
emission  spectra  and  total  primary  photoemlsslon  yields 
for  an  exploding  wire  radiator  (EUR)  source.  The  EUR 
source  of  Interest  Is  generated  by  the  OWL  II*  flash 
X-ray  facility  and  has  been  used  over  the  last  couple 
of  years  by  the  DEFENSE  NUCLEAR  AGENCY  to  Investigate 
SCEMF  response  In  the  soft  X-ray  region. 

The  photoemlsslon  calculations  have  been  carried 
out  using  the  code  SXRP  which  provides  for  a  rigorous 
description  of  electron  transport  to  energies  well 
below  a  kilovolt.  For  a  description  of  the  code,  we 
refer  the  reader  to  the  various  publications  and 
reports  which  document  the  formulation  and  the  various 
applications.1-®  For  the  EWX  source,  we  are  able  to 
compare  our  spectra  with  the  spectrometer  data  of 
Bernstein9  and  our  total  primary  yields  both  with  the 
data  of  Bernstein  and  Fromae,  et  al.1® 

Last  year,  we  reported  on  EWR  photoemlsslon 
results  for  the  materials  At  and  Au. '  At  that  time, 
it  was  not  possible  to  Include  analyses  for  Ag  and  C 
since  the  needed  SXRP  material  models  had  not  yet  been 
developed.  Since  thi*  la  the  fir*t  paper  addressing 
SXRP  results  for  these  newly  modeled  materials,  we 
choose  here  to  also  present  model  information.  A 
fairly  good  description  of  what  a  given  model  contains 
may  be  found  in  last  year's  paper.' 


MATERIAL  MODELS  FOR  Ag  ANO  C 

Given  an  Incident  photon  spectrum,  the  ejected 
electrons  consist  of  two  types  —  photoelectrona  and 
Auger  electrons.  Photoelectrons  from  different  sub¬ 
shells  have  associated  with  them  different  photoionisa¬ 
tion  cross  sections  and  binding  energies.  Subshell 
photolonlsatlon  cross  sections  are  therefore  needed. 
Unfortunately,  compilations  of  theoretical  and  experi¬ 
mental  results  are  all  in  the  form  of  the  total  photo- 
lonlxatlon  cross  section. 11-15  Selected  calculations, 
however,  for  individual  subshells  do  exist.*  “*" 
Complete  compilation  of  the  photolonlsatlon  cross 
section*  for  all  elements,  all  subshells,  and  all 
desirable  energies  are  prohibitive.  We  have  recently 
developed  and  modified  a  code  to  calculate  the  desired 
cross  sections  using  Herman  Sklllman  potentials  and 
wave  functions.  The  performance  of  this  code  is 
satisfactory  In  the  sense  that  the  sias  of  the  subshell 
cross  sections  does  agree  with  experiments  for  photon 
energies  In  the  low  to  sub-keV  range. 

Auger  energies  and  yields  ere  obtained  from 
spectrocopy  experiments1^  and  a  compilation.  Ns  have 
listed  in  Tables  1  and  2  the  binding  energies  and  Auger 
features  for  the  respective  materials,  A-,  and  C. 


For  Inner  shell  Inverse  mean  free  paths  (IMFP), 
we  have  employed  a  technique9  which  provides  a 
relationship  between  the  Impact  IMFP  and  the  corre¬ 
sponding  subshell  photolonlsatlon  cross  section.  In 
the  case  of  conductors,  e.g.,  Ag,  the  IMFP's  for  con¬ 
duction  band  Ionisation  and  plasmon  excitation  are 
obtained  from  the  dielectric  response  function  as 
modeled  by  Ritchie  and  his  colleagues.18* 19  The  IMFP's 
for  Ag  are  shown  in  Figure  1.  The  IMFP's  for  conduc¬ 
tion  band  Ionisation  and  plasmon  excitation  are 
obtained  by  assigning  the  single  0  shell  electron  to 
the  conduction  band.  It  is  Important  that  our  total 
IMFP  curve  agree  with  existing  data,  also  shown  on  the 
same  figure. 

For  carbon,  the  electron  Impact  IMFP's  for  inner 
shells,  a*  calculated  by  the  above  cited  method,  are 
shown  In  Figure  2.  An  Insulator  model  has  been 
developed  by  Ashley,  et  al.19  for  the  IMFP's  of 
polystyrene  and  Its  results  are  also  shown  in  Figure  2. 
The  agreement  Is  satisfactory  considering  the  diversity 
in  the  approaches  used. 

The  elastic  IMFP  is  given  by  the  screened 
Rutherford  formuls.  Values  for  Ag  and  C  may  be  seen 
in  Figure  3  where  they  appear  with  the  total  Inelastic 
IMFP's  from  Figures  1  and  2.  Elastic  scattering  is 
seen  to  be  more  dominant  In  Ag  due  to  Its  higher  2 
value.  For  more  Information  on  our  treatment  of 
elastic  scattering,  we  refer  the  reader  to  the  com¬ 
panion  paper  by  Strickland  and  Lin8  and  to  the  report 
by  Strickland,  et  al.5 


RESULTS 

The  X-ray  source  spectrum  used  In  the  calculation 
of  the  photoemlsslon  yields  is  the  same  one  appearing 
In  last  year's  paper  on  photoemlsslon. i  It  Is  a  con¬ 
tinuum  representation  of  an  actual  OWL  II  X-ray 
spectrum  with  a  prominent  feature  at  M.65  keV  con¬ 
taining  ''.50X  of  the  energy.  For  this  source  represen¬ 
tation,  Figures  *  and  5  show  the  calculated  photoemls¬ 
slon  spectra,  those  measured  by  Bernstein,9  the  calcu¬ 
lated  cumulative  yields,  and  measured  primary  yield* 
by  Bernstein  and  Fromme ,  et  al.10  for  C  and  Ag. 

The  dominant  feature  In  the  C  spectrum  comes 
from  K  photoelectrons  arising  from  Ionisation  by  the 
X-rays  In  the  strong  1.64  keV  feature  in  our  source 
representation.  At  low  energies,  a  feature  appears 
due  to  the  production  of  KLL  Auger  electrons  following 
photolonlsatlon  of  the  E  shell. 

The  dominant  feature  in  the  Ag  spectrum  comes 
from  M  photoelectrons  and  Auger  electron*  following 
photolonlsatlon  of  the  M  and  L  sheila.  The  M  peaks 
correspond  to  the  1.64  keV  source  feature  while  the 
region  above,  until  the  UM  peak  Is  reached,  corre¬ 
sponds  to  the  continuum  region  of  the  source  spectrum 
above  1.6S  keV. 


•This  work  was  supported  by  the  DEFENSE  NUCLEAR  AGENCY 
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The  agreeaent  between  us  and  Bernstein  Is  good 
Cor  C  as  can  be  seen  In  Figure  4.  There  is  also 
reasonable  agreeaent  with  Fronne's  back-biased  diode 
data.  For  Ag,  the  agreeaent  In  spectral  shape  between 
us  and  Bernstein  is  also  good  as  can  be  seen  in 
Figure  S.  Our  yield,  however,  is  about  half  of  both 
Bernstein's  and  Froone's  higher  value.  There  are  a 
nuaber  of  potential  sources  for  this  difference,  such 
as  uncertainties  in  the  X-ray  source  spectrum.  Ue  are 
not  in  a  convenient  position  to  assess  those  uncer¬ 
tainties  associated  with  either  the  source  or  the 
photoemission  experiments.  We  have,  however,  under¬ 
taken  an  assessment  of  the  accuracy  of  our  own  results 
and  this  will  appear  next. 


ACCURACY  ASSESSMENT  OF  THE  Ag  RESULTS 

Unlike  the  previously  Investigated  materials  A £ 
and  Au,  we  have  not  found  published  data  for  Ag  on  the 
primary  photoemission  yield  at  X-ray  energies  falling 
within  the  EUR  energy  region.  We  have  thus  chosen  to 
make  a  consistency  check  between  the  Ag  EWR  results 
and  those  for  Au.  We  have  chosen  Au  because: 

e  we  obtain  good  agreement  with  the 
published  photoemission  data,  and 

e  transport  properties  are  similar 
between  Ag  and  Au,  as  will  be 
demonstrated. 

To  examine  consistency  in  photoemlsslon  from  one 
material  to  the  next,  one  must  compare  photoelectron 
source  spectra  and  transport  properties.  Figure  6 
orovldes  for  a  comparison  of  such  spectra  for  the 
assumed  EWR  source  incident  on  Ag  and  Au.  The  Ag 
spectrum  is  seen  to  be  both  weaker  and  softer.  Its 
Integrated  value  is  1.3  x  104  el/cm3-s  (for  unit  inci¬ 
dent  flux)  and  its  mean  energy  is  1.3  keV  compared  to 
values  of  4.7  x  104  el/cm^-s  and  1.5  keV  for  Au.  If 
transport  properties  were  identical  for  these 
materials,  the  given  electron  source  differences 
would  lead  to  n.3. 3  times  larger  yield  for  Au  compared 
to  the  actual  calculated  value  of  3.6. 

To  examine  differences  in  transport  properties, 
photoemlsslon  spectra  were  obtained  for  both  materials 
with  the  same  electron  source,  namely,  the  Ag  source 
in  Figure  6.  These  results  appear  in  Figure  7  and 
demonstrate  that  transport  characteristics,  for  the 
type  of  source  we  are  considering,  are  quite  similar 
for  Ag  and  Au.  Thus,  differences  in  the  photoelectron 
production  rates  for  these  materials  are  primarily 
responsible  for  differences  in  yields.  From  the 
previous  paragraph,  we  see  that  for  the  assumed 
material  model  parameters,  we  do  have  consistency  in 
photoemlsslon  predictions  between  Ag  and  Au. 

Regarding  model  parameters,  we  believe  our  photolonl- 
zatlon  cross  sections  are  accurate  since  their  sum 
agrees  with  the  published  data.  In  turn,  we  expect 
that  our  calculated  electron  source  spectrum  is 
reasonably  accurate. 


SUMMARY  AND  CONCLUSIONS 

In  this  paper,  we  have  presented  the  material 
model  parameters  for  C  and  Ag  needed  by  code  SXRP  to 
investigate  soft  X-ray  photoealsslon.  Using  these 
models,  photoealsslon  spectra  and  yields  have  been 
obtained  for  an  EWR  source  and  compared  with  the  data 
of  Bernstein  and  Fromme,  et  al.  We  obtain  good  agree¬ 
ment  with  the  spectral  data  for  both  C  and  Ag,  good 
agreement  in  yield  for  C,  but  find  the  yield  data  for 
Ag  to  be  about  twice  the  calculated  value.  A  consis¬ 
tency  test  was  then  applied  to  our  Ag  results  by 
introducing  our  previously  calculated  Au  photoemlsslon 
yield  which  agrees  well  with  the  data.  We  found,  by 
comparing  the  Ag  and  Au  electron  source  spectra  and 
transport  properties,  that  their  yields  are  consistent 
with  one  another  for  the  assumed  material  models.  We 
thus  conclude  that  for  the  assumed  EWR  X-ray  spectrum, 
the  calculated  Ag  yield  should  be  accurate  to  better 
than  the  factor  of  two  discrepancy  exhibited. 

An  estimate  of  the  contamination  effect  for  Ag 
based  on  our  results  of  Al  and  A^O},  indicates  that 
the  photoemlsslon  yield  of  the  contaminated  Ag  will 
be  lower  than  reported  here. 

In  summary,  we  have  compared,  to  date,  our 
photoemlsslon  results  with  the  data  for  an  EWR  source 
incident  on  Ai,  Au,  C,  and  Ag  (Al  and  Au  results 
previously  presented^) .  In  general ,  the  agreeaent  in 
spectral  behavior  between  our  calculations  and 
Bernstein’s  measurements  is  good.  With  the  exception 
of  Ag,  agreement  is  also  satisfactory  between  our 
primary  yield  values  and  the  available  data. 
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TABLE  1.  BIND INS  ENEB6IES  AND  AUGER  FEATURES  OF  Ag 


BINDIN6  ENERGIES  OF  Ag 


SHELLS 

BINDING  ENERGT 
(keV) 

K 

24.? 

L 

3.35 

"l 

0.665 

*23 

0.567 

*45 

0.384 

N 

0.0127 

AUGER  FEATURES  OF  Ag 


ENERGY 

TRANSITION 

(keV) 

TIELD 

KLL 

18.8 

0.166 

LHH 

2.57 

0.900 

NNN 

0.355 

1.000 

TABLE  2.  BINDIN6  ENERGIES  AND  AUGER  FEATURES  OF  C 


BINDING  ENERGIES 


BINDING  ENERGIES 

SHELLS 

(keV) 

K 

0.284 

L 

0.008 

AUGER  FEATURES 


TRANSITION 

ENERGY 

(k«V> 

YIELD 

KLL 

0.272 

1.00 
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FIGURE  S.  CALCULATED  AND  MEASURED 
FMOTO0UMION  SPECTRUM  AND  CUMULATIVE  RACK 
YIELD  FOR  AN  MR  SOURCE  INCIDENT  ON  %. 


FIGURE  7.  PHOTOENISSION  SPECTRA  FRON  Ag  AND  Aw 
FOR  THE  SANE  ELECTRON  SOURCE  SPECTRUM#  NAMELY  THAT 
OF  Ag  SHOMN  IN  FIGURE  6. 
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PhotokMiization  cross  sections,  electron-impact  inverse  mean  free  paths, 
and  stopping  powers  for  each  subshsH  of  silver  a> 

0.  L  Lin  and  0.  J.  Strickland 

Radiation  and  Electromagnetics  Division,  Science  Applications,  Inc,  8330  Old  Courthouse  Road,  Vienna, 

Virginia  22180 

(Received  3  August  1979;  accepted  for  publication  29  October  1979) 

Using  the  Herman-Skillman  potentials  and  bound  wave  functions  for  each  subshell  of  silver,  we 
have  computed  the  continuum  wave  functions,  and  subshell-by-subshell  photoionization  cross 
sections  with  photoelectron  energies  up  to  10  keV.  Applying  a  relationship  between 
photoionization  and  electron  impact  ionization,  we  have  obtained  inverse  mean  free  paths  and 
stopping  powers,  again  by  subshell,  for  electrons  penetrating  through  stiver.  The  maximum 
electron  energy  considered  is  100  keV.  For  the  total  photoionization  cross  section,  comparison  of 
our  work  with  experiment  shows  excellent  agreement  for  photon  energies  down  to  100  eV,  below 
which  solid-state  effects  should  be  included.  Theoretical  total  inverse  mean  free  paths,  being 
strongly  dominated  by  contributions  from  4d  electrons,  are  in  good  agreement  with  data  around  1 
keV,  but  about  a  factor  of  2  larger  at  energies  below  lOOeV.  Our  stopping  power  is  in  good 
agreement  with  other  theoretical  work  above  400  eV  and  approaches  the  relativistic  Bethe 
formula  above  10  keV.  Range  is  also  computed  and  is  in  good  agreement  with  other  theoretical 
work. 

PACS  numbers:  32.80.Fb,  34.80.Dp,  79.20.Kz 


I.  INTRODUCTION 

Interest  in  photon  attenuation  in  materials  dates  back 
to  the  beginning  of  this  century.  The  information  about  at¬ 
tenuation  coefficients  is  frequently  needed  for  the  analysis  of 
radiation  experiments  and  for  application  to  medical,  engi¬ 
neering,  crystallographic,  and  other  practical  problems. 

Most  recently,  soft  x-ray  photoemission  and  electron 
transport  studies  for  radiation  hardening  of  devices'  requires 
detailed  subshell-by  subshell  information  about  photoab¬ 
sorption  coefficients  and  inverse  mean  free  paths  (IMFP’s) 
of  electron-impact  ionization.  Another  area  of  interest  is  sur¬ 
face  physics.  A  technique,  called  electron  spectroscopy  for 
chemical  analysis  (ESCA),  was  developed  in  the  late  1960’s 
in  which  an  x-ray  beam  is  directed  onto  a  solid  specimen  and 
photoelectrons  are  collected  and  their  energies  analyzed. 
Closely  related  to  ESCA  is  Auger  election  spectroscopy 
(AES)  in  which  an  electron  beam  is  incident  on  a  specimen 
and  Auger  electrons  are  collected.  Interpretation  of  these 
spectral  measurements  also  requires  the  above  type  of  de¬ 
tailed  information. 

In  studying  the  penetration  of  charged  particles  into 
materials,  Bethe2  developed  a  simple  stopping  power  formu¬ 
la  applicable  to  high-energy  incident  particles.  In  the  inter¬ 
mediate  energy  region,  the  stopping  power  calculated  from 
the  atomic  picture  is  presumably  adequate.  It  is  generally 
assumed  that  for  electron  energies  below  roughly  100  eV, 
solid-state  effects  will  play  a  dominant  rok  for  stopping.  Re¬ 
cently,  calculations  were  made3  for  the  inverse  mean  free 
paths  and  stopping  power  due  to  conduction  band  ionization 
and  plasmon  excitation  for  various  materials  using  the  meth¬ 
od  of  the  dielectric  response  function.  In  this  paper,  we  try  to 
assess  how  low  in  energy  the  atomic  picture  can  be  applied  in 
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the  calculation  of  subshell-by-subshell  IMFP’s  and  the  stop¬ 
ping  power,  what  the  dominant  contribution  is,  and,  further¬ 
more,  what  kind  of  solid-state  picture  we  should  take  in  im¬ 
proving  our  atomic  calculations  at  low  energies.  In  Sec.  II, 
we  will  briefly  describe  the  theory  of  photoionization  and  a 
relationship  between  the  photoionization  and  electron-im¬ 
pact  ionization.  Section  III  discusses  the  numerical  ap¬ 
proach.  Results  and  comparison  will  be  made  in  Sec.  IV. 

II.  THEORY 

The  transition  rate  for  the  photoionization  process  in 
which  the  incident  photon  is  absorbed  and  an  atomic  elec¬ 
tron  is  ejected  can  be  obtained  from  quantum  theory.  Briefly, 
ad-normalized  continuum  state  with  momentum  k (the  Anal 
state)  can  be  expanded  in  terms  of  spherical  harmonica  and 
radial  wave  functions.  Integrating  over  all  directions  of  k 
and  using  the  asymptotic  behavior  of  the  radial  wave  func¬ 
tion,  one  can  prove  that  the  transition  rate  is  given  by  (wc  use 
atomic  units): 

R  =  lw\<+Ui\H»\il>\\  (I) 

The  continuum  radial  wave  function  Rkl/  of  tfU/  has  the  well- 
known  asymptotic  form 

lr  being  the  angular  momentum  of  the  flaal  electron  and  Hm 
the  interaction  Hamiltonian.  d,  refers  to  the  initial  wave 
function  .Dividing  the  transition  rate  by  theme  id  eat  phoaoa 
flux,  the  photoionization  cross  section  is  obtained: 

|<du,|r«|*<  >  I*-  (2) 

a  is  the  fine-stnicture  constant  and  at  is  the  angular  frequen¬ 
cy  of  the  photon.  The  evaluation  of  the  angular  part  in  tht 
above  matrix  element  can  be  done  by  the  technique  of  re- 
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duced  matrix  elements 
■nation,  the  croas  section  reduces  to 

' —  • 

4 V*a 

& velocity  ““  _  .  ,  ® 

y  3  2/,  +1 

for  length  and  velocity  forms,  respectively.  In  these  expres¬ 
sions,  N„i  refers  to  the  number  of  electrons  in  the  initial 
subshell  nlt .  /„,  is  the  maximum  value  of  /,  and  /,.  P,, 
3irRH,  and  PumtrRUt.  For  radial  wave  functions  in  any  lo¬ 
cal  potential,  these  two  cross  sections  will  yield  identical 
results. 

The  differential  equation  satisfied  by  the  coninuum 
wave  function  Ptt  is  given  by 

(S  +K(r)+  k'~  Pu(r) “  0  (5> 

Once  the  subshell  photoionization  cross  sections  are  known, 
the  corresponding  electron-impact  ionization  cross  sections 
and  hence  IMPFs  can  be  obtained  by  an  approximate  meth¬ 
od  recently  developed  by  us.4  The  formula  we  obtain  is  given 
in  atomic  units  by 

“-Vft— )to— *■  <«> 

iraE,  Jm  \  e  J  q^ 

where  K  is  tk^IMFP’s,E,  represents  the  energy  of  the  inci¬ 
dent  electron,  B  is  the  binding  energy  of  the  subshell  in  ques¬ 
tion,  and  ar(e)  denotes  the  photoabsorption  coefficient  for 
photon  energy,  q^,  and  qmin  are  momentum  transfers: 

*L,  -  IE,  +2<  E,  -  €)  -2  <2£,),'J[2(£,  +  B  -2e)  J*'1, 

(7) 

qL.  =  2 iff,  +2( -  e)  -2  (2£,)l/3[2(£,  -  e)],/J,  (8) 


From  the  subshell  IMPFs,  one  can  compute,  in  atomic 
units,  the  corresponding  stopping  powers 

•m- 

and  ranges 

We  chooae  £q  =  0.37,  corresponding  to  10  eV. 


H(.  NUMERICAL  APPROACH 

To  solve  Eq.  (5),  we  adopt  an  approach  quite  similar  to 
that  of  Cooper5  and  Maneon  and  Cooper.6  The  formulas  for 
the  Runge-Kutta  method  are  obtained  from  Hildebrand,7 
slightly  different  from  Reft.  3  and  6.  The  continuum  wave 
function  Pk,(r)  near  the  origin  is  obtained  from  Hartne.* 
Bound  wave  functions,  central  potentials,  and  binding  ener¬ 
gies  are  tabulated  in  the  literature  *  The  phase  shift  was  com- 
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pitted  following  the  method  suggested  by  Senton  and 

10  ty  nnmnii,«ti«<  co— tint  can  be  calculated  by  * 
««npi>  Mtowinfl  pf  this  method,  as  given  in  the  appendix  rf 
Ref.  3.  The  two  points  chosen  in  the  present  work  for  foe 
calculation  of  foe  normalization  constant  are  usually  sepa¬ 
rated  by  a  few  wavelengths  rather  than  neighboring  grid 
points,  as  adopted  by  others.  To  assure  that  the  correct  nor¬ 
malization  constant  was  obtained,  we  calculated  it  for  a  few 
different  pain  of  points.  Sinoe  we  are  interested  in  cases- 
where  the  energies  of  foe  photoelectron  range  from  a  few  eV 
to  10  keV,  the  iterative  scheme  is  adopted  in  solving  foe 
nonlinear  equation  for  the  amplitude  of  foe  continuum  radi¬ 
al  wave  function,  instead  of  the  first  order  iteration  approxi¬ 
mation  as  used  in  Reft.  3  and  6.  The  iteration  is  stopped 
when  the  absolute  fractional  changes  of  the  quantity  in  ques- 
tion  at  all  points  considered  are  less  than 0.0003.  Usually  two 
iterations  are  needed  for  s$4>  and / photoelectrons  of  silver. 

Calculations  of  the  photoionization  cross  sections  of 
bound  electrons  of  a  given  angular  momentum,  but  different 
principal  quantum  numbers,  were  carried  out  in  one  run.  To 
treat  foe  wide  ranges  in  photoelectron  energies  and  binding 
energies  (hence,  bound  wave  functions),  we  chose  to  have  20 
radial  grid  points  within  a  quarter  wavelength  of  the  photo- 
electron  wave  function  and,  at  the  same  time,  allowed  for  at 
least  200  grid  points  within  the  maximal  radius  of  the  most 
tightly  bound  electron.  Uniform  grid  size  was  doubled  every 
40  grid  points.  The  calculation  of  the  continuum  wave  func¬ 
tion  starts  at  a  radius  where  the  fractional  change 

becomes  smaller  than  one-thousandth  [f  jj(j)  being  the  in¬ 
verse  of  the  amplitude  of  the  radial  wave  function]. 

IV.  RESULTS  AND  DI8CU8SK>N 

The  subshell  photoionization  cross  sections  of  silver 
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Electron  Energy  (eV) 

FIG.  2.  Invene  mean  free  paths  of  electrons  incident  on  Ag.  Remit*  of  two 
other  theoretical  calculation*  (Reft.  1 5  »nd  16)  »re  also  shown  See  Ref.  IS 
for  data  reference* 

calculated  from  the  method  described  above  are  shown  as 
dashed  curves  in  Pig.  1.  The  solid  curve  represents  the  total 
cross  section.  The  data  points  are  taken  from  the  data  compi¬ 
lation  of  Hubbell." 

The  transition  matrix  elements  of  Eqs.  (3)  and  (4) 
change  sign  for  the  case  of  5s  electron  when  the  photon  ener¬ 
gy  passes  through  —  1 1  eV.  This  is  the  so-called  “Cooper 


Electron  Energy  (ev) 

FIG.  4.  Range  of  electrons  in  Aj.  Minimum  electron  energy  »  lOeV. 

minimum” 12  where  the  photoabsorption  coefficient  of  a  par¬ 
ticular  bound  electron  becomes  zero.  The  structure  of  the 
total  cross  section  curve  between  100  and  300  eV  is  due 
mainly  to  the  contributions  from  4 p  and  Ad  electrons.  The  5s 
electron  does  not  have  any  significant  photoionization  cross 
section  in  the  energy  range  shown. 

The  agreement  between  our  atomic  calculation  and  the 
experimental  data  of  solid  silver  is  excellent  in  the  soft  x-ray 
range  of  100  eV  to  30  keV.  Below  100  eV,  we  do  not  antici¬ 
pate  good  qualitative  agreement  between  our  calculations 
and  experiments  due  to  the  neglect  of  the  solid-state  effects, 
which  include  the  formation  of  the  band  structure  and  shift¬ 
ing  of  the  atomic  energy  levels,13  and  the  correlation  effects 
between  electrons.  Qualitatively,  Fig.  1  indicates  the  contri¬ 
bution  at  low  energies  comes  mainly  from  Ad  electrons.  For 
example,  at  the  photon  energy  ~  40  eV,  the  transition  matrix 
element  of  the  Ad  electron  is  about  two  orders  of  magnitude 
larger  than  that  of  5s.  Assuming  no  significant  change  in  the 
subshell  photoionization  cross  section  in  going  from  an 
atomic  picture  to  a  solid-state  picture  the  Ad  band  will  be 
responsible  for  most  of  the  observed  photoabsorption  coeffi¬ 
cient  at  energies  below  iOOeV. 

We  have  also  calculated  the  subshell  photoionization 
cross  sections  using  a  screened  hydrogenic  model. 14  This 
model  is  reliable  for  the  most  tightly  bound  -shell  elec¬ 
trons,  but  gets  worse  as  the  principal  quantum  number  in¬ 
creases.  For  example,  using  a  screening  parameter  of  10.13 
for  the  3s  electron  of  silver,  we  found  the  photoionization 
cross  section  to  be  about  a  factor  of  2  larger  than  that  of  the 
present  approach. 

Relativistic  effects  neglected  in  the  present  calculation 
of  photoionization  cross  section  ate  not  important,  as  isevi- 


dent  when  comparison  b  nude  with  Scofield's  relativistic 
calculation1’  of  Hart ree-Sla ter  subshell  photoionization 
cross  sections  at  1254  and  1487  eV.  Agreement  is  better  than 
5%  for  any  subshell  of  silver. 

Figure  2  shows  the  subshell  inverse  mean  free  paths  for 
incident  electrons  with  energies  ~  10  eV  to  100  keV.  Total 
IMFP’s  and  data  points  are  also  shown  in  Fig.  2.  Our  total 
1MFP  is  seen  to  be  about  a  factor  of  2  larger  than  the  data 
above  a  few  hundred  eV  and  a  factor  of  3  larger  at  lower 
energies.  As  far  as  individual  subshell  contributions  are  con* 
cerned,  the  4 d  subshell  dominates.  Sinee  we  do  not  expect 
good  agreement  between  the  atomic  photoionization  curve 
and  solid-state  photoabsorption  data  at  lower  energies  (tens 
of  eV),  the  results  computed  from  Eq.  (6)  using  this  atomic 
curve  will  not  correspond  exactly  to  the  solid-state  IMFP 
data.  Judging  from  this  consideration,  the  fair  agreement 
mentioned  above  b  viewed  as  satisfactory. 

The  results  of  a  statistical  model,  computed  by  Ashley 
et  al. 16  and  shown  in  Fig.  2,  agree  with  the  data  at  tens  of  eV. 
At  hundreds  of  eV,  our  curve  is  on  the  upper  side,  whereas 
that  of  the  statistical  model  is  on  the  lower  side  of  the  data. 
From  1  to  10  keV,  these  two  theoretical  curves,  based  upon 
quite  different  approaches,  do  agree  within  a  factor  of  2. 

Penn17  has  carried  out  calculations  of  the  total  IMFP, 
again  using  the  method  of  the  dielectric  response  function, 
for  elements  up  to  Z  =  83.  The  energy  range  is  between  200 
and  2400  eV.  For  silver,  his  results  are  also  plotted  in  Fig.  2. 
Being  based  on  the  same  approach,  this  curve  is  quite  similar 
to  that  of  Ashley  el  al. 

Our  subshell  IMFP  curves  can  be  employed  in  Eq.  (9) 
to  generate  subshell  stopping  powers  for  incident  electrons 
with  energies  from  10  eV  to  100  keV.  Figure  3  shows  our 
results.  The  4 d  subshell  dominates  at  energies  below  1  keV. 
However,  the  3 d  subshell  becomes  comparably  important 
for  higher  energies,  due  to  the  large  binding  energy  of  the 
bound  3 d  electrons. 

The  stopping  power  computable  from  the  statistical 
model  is  in  good  agreement  with  our  total  result,  as  shown  in 
Fig.  3.  In  fact,  they  are  within  about  20%  of  each  other  for 
incident  electron  energies  from  1  to  10  keV.  We  have  also 
included  the  stopping  power  computed  from  the  relativistic 
Bethe  formula"1: 

dE  farNe*  mv>E  \ 

~T,  ~  ~^rzr«lW-7V 

-[2(1— 0J),/J  +  l+0J]log2 

+  l-/?2  +  J[l-<l-02)''2l2],  (11) 

where  I,  a  weighted  average  of  the  excitation  and  ionization 
potentials  of  the  atom  of  the  stopping  material,  is  taken  from 
the  semiempirical  formula19 

7  =  9.1Z(I  +1.9Z'2'3)  eV.  (12) 

Our  total  stopping  power  is  in  good  agreement  with  the 
Bethe  formula  at  high  energies  for  which  the  latter  b  applica¬ 
ble.  The  cross-over  energy  of  these  two  stopping  curves  b  at 
~50keV. 

The  range,  as  defined  by  Eq.  (10),  was  also  calculated 
and  is  shown  in  Fig.  4.  For  comparison,  we  also  include  the 
range  curve  from  the  statbtical  model. 16  The  agreement  of 
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these  curves  b  reasonably  good  between  100  eV  and  10  keV 
with  the  crossover  energy  at  ~ 400  eV.  The  relatively  sharp 
rise  of  our  range  near  10  eV  b  due  to  the  small  stopping 
power  there. 

V.  CONCLUSION 

Using  an  atomic  description,  we  have  computed  ab  ini¬ 
tio,  the  subshell  photoionization  cross  sections  followed  by 
inverse  mean  free  paths  and  stopping  powers  of  silver.  A 
range  calculation  was  also  carried  out  u^ing  our  results  for 
the  total  stopping  power.  Theoretical  curves  for  total  quanti¬ 
ties  were  compared  with  available  solid-state  experimental 
data,  as  well  as  other  theoretical  results.  Generally  good 
agreement  b  found  for  energies  (of  either  photon  or  incident 
electron)  above  100  eV.  Below  this  energy,  where  solid-state 
effects  as  well  as  correlation  effects  become  important,  the 
agreement  is  still  fair.  By  separating  the  total  curve  of  a  given 
observable  quantity  into  contributions  from  individual  sub¬ 
shell,  we  are  able  to  assess  the  relative  importance  of  each 
subshell  in  treating  solid-state  effects  or  correlation  effects. 
For  example,  we  find  that  the  4  d  band  of  silver  should  be 
modeled  in  greater  detail  in  the  solid  state. 

Even  though  other  theoretical  works  calculating  sub 
shell  contributions  to  a  given  observable  are  not  available20 
for  comparison,  we  believe  that  our  results  for  individual 
inner  shell  contributions  are  quite  accurate  based  on  the  ob¬ 
servation  that  our  photoionization  cross  sections  are  in  ex¬ 
cellent  agreement  with  data  for  the  inner  shells. 
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